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Abstract 

A thin metamaterial slab comprising a dielectric spacer sandwiched between a metallic grating 
and a ground plane is shown to possess spatially coherent surface resonance states that span a large 
frequency range and can be tuned by structural and material parameters. They give rise to nearly 
perfect angle-selective absorption and thus exhibit directional thermal emissivity. Direct numerical 
simulations show that the metamaterial slab supports spatially coherent thermal emission in a wide 
frequency range that is robust against structural disorder. 
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I. INTRODUCTION 



Surface plasmon polaritons (SPPs) can modulate light waves at the metal-dielectric in- 
terface with wavelength much smaller than that in free space [1.1, which enables the control 

n 

of light in a subwavelength scale for nanophotonic devices [2|j. SPPs with large coherent 
length are useful in many areas, including optical processing, quantum information js| and 
novel light-matter interactions (4|. The enhancement of local fields by SPPs is particu- 



larly important as it opens a new route to absorption enhancement |5y , nonlinear optical 
amplification js, 7] as well as weak signal probingjs, 9[. As the properties of SPP are pretty 
much determined by the natural (plasmon) resonance frequency, there is not much room 
for us to adjust the SPP response for practical applications. With induced surface current 
oscillations on an array of metallic building blocks 10Ml6| , metamaterial surfaces can ma- 
nipulate electromagnetic waves in a similar way as SPPs. Such SPPs or surface resonance 
states on structured metallic surfaces are tunable by geometric parameters. 

In this paper, we examine the properties of surface resonance states at a dielectric- 
metamaterial interface that exhibit magnetic response to the incident waves and strong 
local field enhancement. We will see that these surface resonance states can give highly 
directional absorptivity and emissivity, and may thus help to realize interesting effects such 
as spatially coherent thermal emission. As the structure is very simple, it can be fabricated 



down to the IR and optical regime 17H19|. 



We will show that a thin metamaterial slab, with a thickness much smaller than the 
operational wavelength, supports delocalized magnetic surface resonance states with a long 
coherent length in a wide range of frequencies. Operating in a broad frequency range, 
these spatially coherent SPPs are surface resonance states with quasi-TEM modes guided 
in the dielectric layer that are weakly coupled to free space, and the coupling strength 
can be controlled by tuning structural parameters while the frequency can be controlled 
by varying structural and material parameters. The high fidelity of these surface resonance 
states results in directional absorptivity or emissivity, which is angle- dependent with respect 
to frequency. Finite-difference-in-time-domain (FDTD) simulations verify that the highly 
directional emissivity from the slab persists in the presence of structural disorder in the 
grating layer. 

Such metal-dielectric-metal (MDM) structures were recognized as artificial magnetic sur- 
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12|]., the magnetic response were 



20| . After the concept of meta- 



faces with high impedance by the end of last century 
described with an effective permeability in Lorentz typef 

material being proposed}^. P. Alastair and his co-workers numerically and experimentally 
proved that the ultra-thin MDM structures can resonantly absorb or transmit radiations at 
low frequency limit [22]. They addressed that the central frequencies of absorption peaks are 
independent from the incident angle with an interpretation of Farby-Perrot resonant mode 
(EQ. 1 in Ref. 22). The same group further explored the angle-independent absorption, 
as the main scenario of the incremental work, by measuring the flat bands of surface wave 



dispersion in the visible [231] as well as the microwave region [24]. In contrast, we find that the 
structure with proper design also supports very narrow absorption peaks which are sensitive 
to the incident angle and obviously do not satisfy to the Fabry-Perot resonance condition 
suggested in the previous studies. 

It is worth noting that the physics origin of an angle-independent peak is quite different 
from that of angle-dependent ones. The former, investigated in Refs. 22-24, mainly comes 
from the localized surface resonance states, while the latter, found by us, comes from collec- 
tive surface resonance states. An intuitive picture is as follows: high order quasi-TEM modes 
induced inside the dielectric layer can assign phase correlation to the outgoing waves emitted 
from the air slits of grating, thus are very crucial to the formation of collective response. 
Weak enough both the leakage from dielectric layer to air slits and the material absorption, 
the spatial coherence of surface resonance states will survive. As the interaction between 
the structure and the incident waves will excite quasi-TEM modes inside the dielectric layer, 
the magnetic induction must be parallel to the MDM surfaces if it exists. Thus a surface 
resonance state on a MDM structure is usually magnetic in nature. Our findings about 
spatially coherent surface resonance states are original compared to the common knowledge, 
and have great potentials in coherent control of SPPs as well as thermal emission radiations. 



II. MODEL DESCRIPTION AND MODE EXPANSION METHOD 

Our model system is schematically illustrated in Fig. 1. Lying on the xy plane, the slab 
comprises an upper layer of a metallic lamellar grating with thickness t, a dielectric spacer 
layer as a slab waveguide with thickness h and a metallic ground plane. The metallic strips 
are separated by a small air gap g, giving rise to a period of p = a+g for the lamellar grating. 
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FIG. 1: Schematic picture of the magnetic metamaterial slab. The geometric parameters are 
t = 0.2fim, h = 0.8/xm, a = 3.8pm, g = 0.2//m and t p = a + g = 4.0/xm. the dielectric spacer 
layer is slightly dissipative by assigning a complex permittivityej-j-j- = e r e^ + i^- with e r = 2.2 
ander = 66.93S/m. 



Each metallic strip together with the ground plane beneath it constitutes a planar resonant 
cavity as the building block that gives magnetic responses at cavity resonances 12|, |l6) . The 
metallic grating is along the x direction so that the guided waves in the dielectric layer (at 
< z < h in region III) can only couple with the transverse magnetic (TM) polarized wave 
(with the electric field-Ein xz plane) in the free semi-space (at z > h + 1 in region I). The 
geometric parameters of our model are t = 0.2/im, h = 0.8/im, a = 3.8/im, g = 0.2/im 
and p = a + g = 4.0/im. For a TM incident plane wave with an in-plane wave vector 
kon = k x e x + k y e y , the EM field in region I and in region III can be written in terms of the 



reflection coefficients r m and the guided Bloch wave coefficients t™ |25 



as 
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where the term 5 m o e 



-ikl 



(r\ko..) denotes the incident plane wave with (rl&o,,) 



e i{k x x+kyy) ^ m0 being the Kronecker function and m being the Bloch order; (rj^) 



MO 



denotes wave component of the m Bloch eigenmode in the semi-free space 



(region I) and the dielectric layer (region III) with respect to k m = ko« +G m . k m is the in-plane 
wave vector and G m = ^^e^ is the m th reciprocal lattice vector. k J z = y EqIIqUJ 2 — \k m \ 2 



and k 1 ^ = y eni^o^ 2 — \k m \ 2 are the z components of wave vector for the m th order Bloch 
eigenmode in region I and region III respectively. Eq and Em are the permittivity of the 
vacuum and the dielectric slab, /i is the vacuum permeability. 

We shall mainly consider infrared frequencies, at which the metals can be well approxi- 
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mated as perfectly electric conductors (PEC). The EM fields a.th<z<h + tm region II 
are squeezed inside the air gaps, in which the magnetic fields can be expressed in terms of 
the expansion coefficients a; and b\ of forward and backward guided waves, as: 

H 2 (f, z) = J2 W~ iqi{z ~ h ~ t) + b ie iqi{z - h) } (r\ai) , (2) 
i 

where {r\ai) = cos[lir/g(x + g/2)], (I = 0, 1, n, ...) is the in-plane distribution of guided 
mode \oti) running over all air gaps defined by 29]. qi = yj 'e fi Q u 2 — (lir/g) 2 — k 2 is the z 
component of wave vector for the I th guided mode \oti). 

We can obtain the coefficients t m (f, k 0// ) and r m (f, k 0// ) of the m th guided and reflected 
waves by applying the boundary continuity conditions for the tangential components of 
electromagnetic wave fields (over the slits) at the interfaces z = h and z = h + t. Given 
that surface resonance modes are intrinsic response, we can also assign zero to the incident 
plane wave and apply the boundary continuity conditions for the tangential components of 
wave fields to derive the eigen- value equations. A surface resonance state can be determined 
by searching a zero value/minimum of eigen-equation determinant in the reciprocal space 
provided that it is non-radiative/radiative with infinite/finite life time below/above light 
line in free space. 



III. ABSORPTION SPECTRA PROPERTIES AND SPACIAL COHERENCE OF 
MAGNETIC SURFACE RESONACE STATES 

We derived the absorption spectra of the slabv4(fc , = 1 — -^(tt 2 ") \ r m(ko, v) | 2 which 

m u 

includes the contributions from all Bloch orders of reflected waves. As a consequence, 

— * 

A(ko, uj) gives information about the surface resonance states as well as the emissivity prop- 
erties as governed by Kirchhoff's lawl^] . We shall assume that the dielectric spacer layer 
is slightly dissipative by assigning a complex permittivity^/// = e r Eo + ia/u with e r = 2.2 
and a = 66.93S/m [Im(ejjj) ~ 10~ 2 e r e ] in the calculated frequency regime. In Fig. 2(a), 
we present the absorption spectra at various incident angles. The spectra exhibit a low 
and broad peak at 13.2THz which is almost independent of the incident angle, while the 
other absorption peaks at higher frequencies are narrow and sensitive to the incident angle 
with a maximum absorption approaching 100%. The slab thus acts as an all-angle absorber 
at 13.2THz (a similar result can be found in Ref. 31), but exhibits sharp angle-selective 
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FIG. 2: Absorption spectra under TM-polarized incidence (a) as a function of frequency at incident 
angles of 9 = 0°, 5°, 20°, 30° (g=0.2//m) and (b) as a function of incident angle at 54.3THz (solid 
line, g=0.2fim), AOTHz (dashed line, g=0.2^ra) and 40THz (dotted line, g=0.1/um) 

absorption peaks at higher frequencies. Shown as solid and dashed lines in Fig. 2(b), the 
sharp angular dependence of absorption coefficients (note that the vertical axis is in log- 
scale) at 40.0THz and 54.3THz implicitly implies the existence of spatially coherent surface 
resonance states. The angle- dependent absorption peaks become lower and disappear grad- 
ually with the increase of the material loss. This presents a way to realize nearly perfect 
absorption with weakly absorptive materials by coherent surface resonance states. The co- 
herent length of a surface resonance state can be estimated by the ratio of the wavelength 
A and the full width at half maximum (FWHM) A9 of the absorption peak 32| . For ex- 
ample, for the T 4 state at 54.3THz and ko.. = 0, the angular FWHM of the corresponding 
absorption peak AO = 4.6°(from 9 = — 2.3°to 9 = 2.3°) gives rise to a coherent length 
X/A9 = 68.5/im « 12. 4A. The coherent length is about 220A for the surface resonance state 
at 50.22THz and k « 0.02vr/p with A9 = 0.26°(not shown in figure). The angular FWHM 
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FIG. 3: (a) Dispersion diagram of TM-polarized surface resonance states, (b) The reflection phase 
difference between the th order reflection and the TM-polarized incident plane wave at incident 
angles of = 0°(red line) and 9 = 2° (blue line). 

is reduced if the gap size is smaller, as shown with the dashed and dotted lines in Fig. 2(b) 
for g = 0.2/im and g = 0.1/xm at 40THz, which means that the coherent length of the surface 
resonant modes can be controlled by the gap-period ratio g/p. 

To quantitatively characterize the formation of these spatially coherent surface resonance 
states, we employ the eigenmode expansion method to calculate the surface resonance disper- 
sion (in the limit of no material loss) as shown in Fig. 3(a). The B\ surface resonance states 
lie blow the light line L 2 (magenta dashed line), and thus are non-radiative as evanescent 
modes. The surface resonances labeled as B 2 originate from the coupling of the fundamen- 
tal magnetic resonance modes of the metal strip structure with the free space light line L 2 . 
The surface resonances B 3 and B± are harmonic modes of the magnetic resonances that 
hybridizes with the guided mode inside the dielectric layer. The calculated reflection phase 
difference between the th order reflected and incident electric field, as shown in Fig. 3(b) 
for normal incidence (red line), and 2° incidence (blue line), clearly shows that the resonances 
are magnetic in nature when the surface resonances intersect the zone center at T2(13.2THz 
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FIG. 4: Spatial distributions of magnetic fields and electric fields in the xz plane for T 2 state at 
fr 2 = 13.2THz, 6 = 0°[(a) and (c)], a state on B 3 at / = 50.22THz,<9 = 2° [(b) and (e)]and T 4 state 
at fr 4 = 54.3THz, 6 = 0°[(c) and (f)] 

) and r 4 (54.3THz) as the reflection phase is zero like what a magnetic conductor surface 
does to the incident waves. The state r 3 , invisible in the reflection phase spectrum under 
normal incidence [red solid line in Fig. 3(b)], is a dark state as its eigenmode is in mirror 
symmetry about the yz plane and can not couple with free space photons. While the other 
B 3 states can couple with external light under oblique incidence [see the blue line in Fig. 
3(b)]. For example, there exists in-phase reflection at frequency 50.22THz under an incident 
angle of 2°, corresponding a B 3 state at frequency 50.22THz and fco,, ~ 0.027r/p. 

The angle-independent absorption peak at 13.2THz is due to the B 2 mode, which is 
only weakly dispersive near the zone center. The more dispersive B 3 and B^ modes are 
accountable for the incident-angle sensitive absorption in the higher frequencies in Fig. 
2(a). The field patterns in Figs. 4(a)-4(c) present the spatial distribution of the real part 
of magnetic fields excited by the incident plane waves with incident angles 0°,2° and 0° for 
the three surface resonance states on B 2 , B 3 and B4 respectively, and the corresponding 
vector diagrams of electric fields are shown in Figs. 4(d)-4(f). We can see clearly that the 
electric fields reach maximum in strength at the slab upper surface, and exponentially decay 
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FIG. 5: Magnetic field H inside the dielectric layer normalized to that of incidence Hq under an 
incident angle of 9 = 2°. Black line: all Bloch orders of TEM guided modes included; Red line: 
only the 0*^ Bloch order considered; Blue line: summation of -1 st and +l st Bloch orders of TEM 
guided modes. 

along the surface normal into the free space. This is precisely a picture of SPP modes. The 
field patterns comes from the coincidence of the evanescent wave components in high Bloch 
orders at both sides of metallic grating. 

Although only TEM guided modes are allowed to be excited in the thin MDM slab 
within the frequency of our interest, the B 2 states are quite different from the B 3 and B4 
states in field patterns inside the dielectric layer. We see from Figs. 4(b), 4(c) ,4(e)and 
4(f) that for a B 3 or _B 4 state, there are nodes and anti-nodes in field patterns, while for 
the B2 state, the magnetic field is almost uniformly distributed. Calculations on local field 
enhancement inside the dielectric slab resolve the puzzle. Black solid line in Fig. 5 presents 
the normalized magnetic field \H\ inside the dielectric with respect to that of incidence \H \ 
under an incident angle of 6 = 2°. The Fourier component in m = order [blue solid line] 
contributes the most at 13.2THz and the least at 50.22THz and 54.3THz; while it is just 
opposite for the contributions in combination from the two high order Fourier components 
with m—zkl ordersfred solid line]. Figure 5 also indicates that the enhancement of local field 



9 



of an excited B 3 or B A state can be ten times larger than that of an excited B 2 state; the 
enhancement factor at 50.22THz is about 100, while it is only 10 at 13.2THz. 

We see from Fig. 3(a) that the surface resonance dispersion of the slab comes from the 
interaction between the magnetic resonances and the (folded) light lines L\ (for dielectrics) 
and L 2 (for air) grazing on the interfaces. In the limit of a small gap-period ratio (g/p = 0.05 
for example), our system is weakly Bragg-scattered, and as such, when a surface resonance 
state on branches B 3 or B 4 is excited, the induced wave fields inside the dielectric of region 
III are guided quasi-TEM modes dominated by ±l st Bloch orders. For that reason, the 
i? 3 and B 4 states have high fidelity even though they are leaky modes, as most of their 
Bloch wavefunction components lying outside the free space light line. As the air gaps of 
the metallic grating serve to couple the electromagnetic waves of region I and region III, 
the quality factor of a resonance state can be estimated with the overlap integral between 

— * . 

the fundamental waveguide mode \a ) in the air gap and the dominant Bloch waves \k l m ) 
[i = I, III) in region I or region III for the coupling coefficients 

C' m = ( Q „|fl^> = -i,^^^ sinc[*^=^], (3) 

when the air gap width g <C p is satisfied. For the B 2 states, the major Fourier component 
of the wavefunction is \k l ) m zero order, and as k 1 ^ 1 is generally not small, Cq 11 is usually 
very large according to Eq. 3, and the B 2 states leak out easily. The states on branches B 3 
and B 4 have major Fourier components in m=±l order, and as they are asymptotic to the 
(folded) dielectric light lines Li, the absolute value of k 1 ^ (m = +1 for k x < or m = — 1 
for k x > 0) is very small, resulting in the small coupling coefficients C u l or C+y 1 . The B 3 
and _E>4 modes have to travel a long distance before they leak out. They have a long life time 
and a good spatial coherence. It also explains why the state T 3 , a state precisely superposing 
on folded light line L[ in dielectric layer, is dark to the incident plane wave as = 0. 

Different from B 3 and _B 4 states, the B 2 states have a major Fourier component in m = 
order which directly couples to the free space photons. As a consequence, the B 2 states, 
forming a flat band far away from the light line L 2 when fc ,, is small, are localized with 
resonant frequency scaled by local geometry of unit cell. The high model fidelity of a B 3 or 
-B4 state also gives rise to much more intense local field compared to the B 2 states. As shown 
in Fig. 5, the induced local field is 100 times stronger than the incident field for the state 
on B 3 ; while it is only 10 times stronger for r 2 , and this is consistent with the absorptivity 
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shown in Fig. 2(a). In addition, the coherent length can be adjusted by the gap width as the 
kernel (ao\r\k l ) is proportional to the gap-period ratio g/p. More calculations demonstrate 
that the angular FWHM of the absorption peak is reduced from 0.26° to 0.16° when the gap 
is decreased from 0.2/zm to 0.1/im, corresponding to a coherent length of 358A. 

We note that most of the attentions in previous studies have been devoted on the localized 
22IJ24, 



B 2 states 



31] . While the spatially coherent surface resonance states will lead us into 



a new vision about coherent control of emission radiations. J. -J. Greffet and co-workers 
showed that highly directional and spatially coherent thermal emission can be obtained by 



32M35| . The magnetic resonant modes 



etching a periodic grating structure into a SiC surface 
in our system can do the same, as will be demonstrated below. Our system has the advantage 
that the operational frequency is tunable by changing the structural parameters, and the 
operational bandwidth is wide. In addition, our structure supports all-angle functionality 
for some specific range of frequencies as shown in Fig. 6, although it is periodic only in one 
direction. 



IV. COHERENT THERMAL EMISSION 

We performed finite-difference-in-time-domain (FDTD) simulations to emulate the emis- 
sions from a slab containing point sources with random phases using the same configuration 
parameters aforementioned. We purposely put disorder in structure to test the robustness of 
the phenomena. We assigned two Gaussian distributions(they can be uniform distributions 
or other types as well) independently to the width of metallic strips and the center posi- 
tions of air gaps to introduce a 4% (standard deviation) structural disorder. The slab has a 
lateral size of 60 periods along the x direction. A total of 1200 point sources with random 
phases are placed at the mesh points inside the dielectric layer. Directional emissions of a 
wide range of frequencies above 34THz are confirmed by the simulation. The 4% structural 
disorder has little impact on the directional emissivity. Figs. 6 (a), 6(c) and 6(e) show the 
far-field emission patterns in the xz plane (H-plane) at 40.0THz, 54.3THz and 58.0THz. 
The inset in Fig. 6(c) is a control calculation in which the top metal gratings are removed, 
so that there is just a dielectric layer with random phase sources above a metal ground 
plane. The directivity of emission from the random sources is lost. Figs. 6(b), 6(d) and 6(f) 
present the absorptivity (under plane wave incidence) as a function of in-plane wave-vector 
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FIG. 6: Radiation patterns in the H-plane (calculated by FDTD) and absorptivity (calculated 
by mode expansion method) as a function of in-plane wavevectors at / = 40.0THz[(a) and (b)], 
/ = 54.3THz[(c) and (d)] and / = 58.0THz[(e) and (f)]. In FDTD simulations, 1200 point sources 
with random phases are placed at the mesh points inside the dielectric layer. A 4% structural 
disorder is included in the80/misimulation cell, which accounts for the slight asymmetry of the 
radiation patterns, but also demonstrates the robustness of the angle selectivity with respect to 
disorder. The inset in (c) is a control calculation in which the top metal gratings are removed, so 
that there is just a dielectric layer with random phase sources above a metal ground plane. The 
directivity of emission from the random sources is lost. 
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ize coherent thermal radiations, such as 
36| or one-dimensional photonic crystal 



(solid angle) at these frequencies. The strong angle selectivity of the absorption is evident, 
and by Kirchhoff's law, the thermal emission should also be highly directional, which is a 
direct consequence of the good spatial coherence of the surface resonance states. As shown 
in Figs. 6(b), 6(d) and 6(e), the absorption/emission peaks generally trace out an arc in 
the k x ~ k y plane, but near 54.3THz[Fig. 6(d)], the dominant emission beam is restricted 
to a small region near the zone center. This is because the 1?4 state is a minimum point if 
we consider the band structure in the k x — k y plane. That means that at 54.3THz, we can 
obtain a directional emission beam not just in the H-plane, but in all directions, although 
the structure is periodic in only one direction. 

We note that there are other schemes to rea 
by utilizing three-dimensional photonic crystals 
cavities 37J • Our metamaterial slab presents a route to achieve linearly polarized coher- 
ent thermal emission radiations in a wide frequency range which can be tuned by adjusting 
structural parameters and material parameters. 

V. CONCLUSION 

In summary, we proposed a simple metamaterial slab structure that possesses spatially 
coherent magnetic surface resonance states in a broad range of frequencies. These states 
facilitate nearly perfect absorption in a thin metamaterial slab containing slightly absorptive 
materials. As the absorption spectrum is highly angle- selective, the slab should give direc- 
tional thermal emission. Direct FDTD simulation with random-phase sources corroborates 
the existence of strong angular emissivity even in the presence of structural disorder. As 
the surface resonances originate from artificial resonators, the operational frequency and the 
response can be tuned by varying the structural configurations. The simple metamaterial 
structure may be a useful platform to realize the coherent control of thermal emissions, 
optical antennas, infrared or THz spectroscopy as well as photon detector. 
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